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Abstract. New data from HERA experiment on deep inelastic scattering have 
been used to parametrize nucleon and Pomeron structure functions. Within 
the Gribov theory, the parameterizations were employed to calculate gluon 
shadowing for various heavy ions. The latter was compared with predictions 
from other models. Calculations of multiplicity reduction due to gluon shad- 
owing for d+Au collisions at forward rapidities at ^/snn— 200 GeV are in good 
agreement with BRAHMS data on the nuclear modification factor. 
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1. Introduction 

For low energies of the incoming beam in a hadron-nucleus collision, the successive 
elastic rescatterings of the initial hadron on the various nuclei of the nucleus are 
well described within the probabilistic Glauber model At higher energies, 

corresponding to E crit ~ mjyfiRA, the hadronic fluctuation length can become of 
the order of nuclear radius and there will be coherent interaction of constituents 
of the hadron with several nucleons of the nucleus. Within the Gribov approach [ 
E], this corresponds to summing up contributions of inelastic intermediate states, 
and leads to a reduction of the total cross section of the reaction, i. e. to nuclear 
shadowing. 

We calculate the total amount of gluon shadowing for low values of the Bjorkcn 
variable x for heavy ions, ignoring for the time being the contribution from the 
quarks. The most recent data on diffractive structure functions are used and much 
stronger shadowing effects than previously expected are found. These effects will 
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lead to a strong multiplicity reduction in A+A collisions at RHIC and LHC energies. 
2. The Model 

The diffractive j*N scatterings are described by Pomeron exchange. The scattering 
amplitude of an incoming photon with virtuality Q 2 on a nuclear target, consisting 
of A nucleons, can then be written as 

a A = Aa N + cr^ 2) + ... . (1) 

The second term in £Q) is negative and is related to diffractive DIS through the 
AGK cutting rules Higher order rescatterings in are model dependent. 

The Schwimmer unitarization [0 for the 7* A cross section is used to obtain the 
shadowing ratio 

R Sch (x) = a i*A _ f d 2r T A{b) , - 

where f{x,Q 2 ) is the effective shadowing function, Ta(6) is the nuclear density 
profile normalized to unity and standard DIS variables are used. Following [0 [7] 
in choice of parameters and kinematics, one can get the shadowing function as 

f(x,Q 2 ) = ^ r dxjpBixjp) ^T'^:® Flitmin.) ■ (3) 



F2(x,Q 2 ) 

Here B(xp) = 0.184— 0.02 In (xjp) fm 2 , and Fa is the form factor of the nucleus. 
Calculations are made both for Xp ax = 0.1 as in [ |S] and for Xp ax = 0.03 

as in [|Hj. The structure functions F2(x,Q 2 ) and F^j {xp , Q 2 , (3) are determined 
from experiment. At small x, gluon shadowing is found to be dominant. Quark 
contribution to the structure functions is not considered in what follows. Shadowing 
due to quarks, obtained wihtin the same approach, was discussed in 

The gluon parton distribution functions (PDF) for nucleon and Pomeron were 
measured at intermediate Q 2 at the HERA experiments ZEUS and HI, correspond- 
ingly. The next to leading order (NLO) ZEUS-S QCD fit for the gluon PDF of 
the nucleon at Q 2 = 7 GeV 2 , and the gluon PDF for the Pomeron (diffractive 
structure function) [0 at Q 2 = 6.5 GeV 2 were both parametrized by 

xg(x,Q 2 ) = Ax- 5 {l-xf , (4) 

where the fitting parameters {A, S, 7} = {1.9,0.19,6.7} were obtained for the nu- 
cleon and {0.38,0.28,0.17} for the Pomeron case, respectively. The Q 2 -dependence 
of the fitting parameters is weak for moderate Q 2 , and so we neglect it for the sake 
of simplicity. 
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Fig. 1. Gluon shadowing for heavy ions. Fig. 2. Comparison of theoretical pre- 
Closed (open) symbols are for x r p ax = 0.1 dictions for the Pb/nucleon ratio, at fixed 
(0.03). Q 2 . 



3. Numerical results 

Gluon shadowing for various heavy ions (Ca, Pd and Pb) from J5J is presented 
in Fig. n The gluon shadowing is very strong at small x, and disappearing at 
x = Xp ax . This is a consequence of the coherence effect in the form factor, and 
the vanishing integration domain in ©. Gluon shadowing is as low as 0.2 for the 
Pb/nucleon ratio. 



A comparison of our results for Pb/nucleon ratio at Q 2 — 6.5 GeV with : 



0.03, with those of others, calculated at Q 2 — 5 GeV 2 is presented in Fig. [21 For x < 
10~ 3 our model predicts the stronger gluon shadowing compared to [^j (dashed- 
dotted line) and (dotted line), while [^| (dashed line) predicts the strongest 
effect down to x ~ 10 -4 . 



4. Shadowing effects in d+Au collisions 

The model is now employed to study multiplicity reduction in d+Au collisions at 
ultra-relativistic energies. Deuteron is treated as a point-like particle in impact 
parameter space, but with the shadowing found from J5J. The collision is described 
by two-jet kinematics through x p (t) = cpr e ±v / \fs, where pr is the transverse 
momentum of the particle, and fixed at Q 2 . We assume that most of the high- 
er particles come from jets c times more energetic than the measured one. The 
theoretical prediction is given by [ I12| 



it. 



theo 
d+Au 



= R s d ch {x p )Rl 



(•5) 



From here one obtains the multiplicity reduction due to shadowing compared to the 
Glauber model. Then the model predictions for nuclear modification factor (NMF) 
at forward rapidities are compared to BRAHMS data at sJsnn = 200 GeV [E]. 
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Fig. 3. NMF ratio for (a) c = 3 and (b) c = 5. See text for details. 



We exclude the gluon shadowing effects in the BRAHMS NMF at r\ = by 
defining i?™^™ = [-R^u/^d-Au] =0 - The multiplicity reduction due to shadowing 
effects will appear when we compare the NMF at forward rapidities, r\ = 1, 2.2, 3.2 
to The ratio R = \R e JJ u ] / 'R^aT is plotted in Fig. together with the 

predictions of (JSJ) for two different values of the parameter c. Statistical errors are 
denoted by the thick solid line, while the systematic and statistical errors added up 
are denoted by the dashed line. Cronin effect is assumed to be rapidity independent 
and is cancelled out in the ratio. The choice of c does not affect the result. Within 
the presented framework, one can conclude that suppression of the nuclear modifi- 
cation factor at forward rapidities is mostly due to gluon shadowing in the nuclei. 
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